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Dyson's Brownian motion model is a one-parameter family of log-potential interacting 
particle systems in one dimension parametrized by an inverse temperature /3 > 0. When j3 = 
1, 2 and 4, this model can be regarded as a stochastic realization of the eigenvalue statistics 
of Gaussian random matrices. Dunkl processes are mathematically defined using differential- 
difference operators (Dunkl operators) associated with finite abstract vector sets called root 
systems. When the root system is specified to be of type A, Dunkl processes constitute a one- 
parameter family of interacting particles in one dimension, in which particles perform not 
only diffusive motion and mutual repulsion but also interchange positions spontaneously. 
In the present paper, we prove that the type-^ Dunkl processes with parameter fc > 
starting from any symmetric initial configuration are equivalent to Dyson's model with the 
inverse temperature /3 ~ 2k. We focus on the intertwining operators, since they play a 
central role in the mathematical theory of Dunkl operators, but their general closed form 
is not yet known. Using the equivalence between symmetric Dunkl processes and Dyson's 
model, we extract the effect of the A-type intertwining operator on symmetric polynomials 
from these processes' transition probability densities. In the zero-temperature limit, the 
intertwining operator maps all symmetric polynomials onto a function of the sum of their 
variables. This allows for the analysis of the zero-temperature limit of Dyson's model, which 
becomes a deterministic process with a final configuration proportional to a vector of the 
roots of the Hermite polynomials multiplied by the square root of the process time, while 
being independent of the initial configuration. 
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I. 



INTRODUCTION 



Dyson's Brownian motion model was introduced as the stochastic process that the eigenvalues 
of a Gaussian random matrix undergo when its entries are independent Brownian motions [1] . In 
this sense, Dyson's model is a dynamic version of Gaussian random matrices [2], and therefore 
both models share the same parameter /3 for each of the three Gaussian ensembles of random 
matrices, 1, 2 and 4 for the Gaussian orthogonal, unitary and symplectic ensembles respectively 
[2-4]. While /3 is usually considered as a discrete parameter, in the case of Dyson's model it can 
be regarded as a real positive parameter. In that case, Dyson's model is a system of N Brownian 
particles in one dimension that repel each other via a log-potential, and 13/2 is its coupling constant 
[5]. Additionally, in view of the fact that the joint probability density function of the eigenvalues 
of Gaussian random matrices has the form of the partition function of a one-dimensional log-gas 
[3], P also carries the meaning of inverse temperature. The vicious walker model formulated by 
Fisher [6] consists of N random walkers who annihilate each other if they make contact; one of the 
most particular properties of Dyson's model is that for [3 = 2, Dyson's model is a diffusion scaling 
limit of the vicious walker model with the restriction that all walkers must survive at all times [7] . 
This restriction is equivalent to conditioning the particles never to collide, and therefore this 
process is also known as the non-colliding Brownian motion [8, 9]. Due to this connection, Dyson's 
model has found applications in polymer networks [10, 11], traffic models [12], gauge fields and 
nuclear physics [13-17], and random growth models [18, 19], among others. In particular, Dyson's 
model has been applied extensively in the Kardar-Parisi-Zhang universality class [20-23] , and there 
have been recent experimental observations of phenomena pertaining to it that confirm theoretical 
predictions [24, 25]. Also, because of this relationship, Dyson's model is intimately related to the 
mathematical areas of combinatorics, group theory and representation theory [26-29]. Therefore, 
Dyson's model has been an active topic of research since its introduction, and remains one of the 
most important models in non-equilibrium statistical mechanics and in stochastic process theory 
as well as a topic of great interest in mathematics. 

On the other hand, Dunkl operators are differential-difference operators [30] that depend on a 
series of parameters and a finite set of vectors. This set of vectors is called root system. In this 
paper, we will focus on the ^-type root system and its corresponding Dunkl operators given by 



[31] 




N 



(1) 
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where A; > is a real scalar parameter, a; is a vector in M^, aijX denotes the vector x with its i-ih. 
and j-ih. components interchanged and the function / is an arbitrary differentiable scalar function. 
This and other types of Dunkl operators have been used in mathematical physics in the study of 
integrable systems and quantum many-body models [32-34], but the examples of applications of 
these operators are scarce in comparison to Dyson's model. While initially defined as a tool to 
study multivariate orthogonal polynomials [30], Dunkl operators have also been used in harmonic 
analysis [35] and stochastic processes [36]. 

Dunkl processes are defined as follows [37]. Consider first the iV-dimensional heat equation. 
The AT-dimensional heat kernel is the transition probability density (TPD) that corresponds to 
the heat equation as a Kolmogorov backward equation. Therefore, one can use the heat kernel to 
define A/^-dimensional Brownian motion. Now, replacing the spatial partial derivatives in the heat 
equation with Dunkl operators such as the one given in Eq. (1) yields the Dunkl heat equation. 
Then, the TPD that solves the Dunkl heat equation as a Kolmogorov backward equation defines a 
generalization of Brownian motion. Dunkl processes arc defined as this generalization, that is, as 
the family of stochastic processes that obey the TPD that corresponds to the Dunkl heat equation. 
Hence, different root systems define different kinds of Dunkl processes. In the case of the type- 
A Dunkl operators (1), the corresponding ^-type Dunkl processes are a one-parameter family of 
stochastic processes, where the real parameter k is non- negative [37]. Dunkl processes subject to 
symmetric initial conditions with respect to their respective root systems are called radial Dunkl 
processes [38, 39], and their properties have been an active topic of research in recent years [40, 41]. 
In order to emphasize their variable exchange symmetry, we will refer to the ^-type radial Dunkl 
processes as symmetric Dunkl processes. 

Many of the results from the theory of Dunkl operators are achieved through the use of Dunkl's 
intertwining operator [42], which we review in Section III. Consider a function f{x) as in 
Eq. (1). Then, Vk is defined by the following equation: 

TiVkfix) = Vk-^f{x). (2) 

In other words, binds Dunkl operators with regular partial derivatives. Additionally, it has 
no effect on constants (i.e., V^l = 1) and when it is applied on a monomial of a given degree, 
it produces a homogeneous polynomial of the same total degree. Note that there is a different 
operator for each type of Dunkl operator. Through the use of Vk, many calculations are simplified, 
as it allows Dunkl operators to be treated as partial derivatives. However, it presently has no 
general closed form. 
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After reviewing Dyson's model briefly in Sec. II A, we will show in Theorem 1, Sec. II B, that 
the TPD of symmetric Dunkl processes obey a Kolmogorov backward equation that is identical to 
that of Dyson's model. Therefore, their TPDs are equivalent and a relationship between the Dunkl 
parameter k and the parameter /3 in Dyson's model is established. This relationship provides a 
physical meaning to the abstract Dunkl parameter fc, that is, k can be understood as a parameter 
proportional to the inverse temperature. In Sec. Ill, we will review the intertwining operator and 
its defining properties. Using Theorem 1, we will extract the effect of the intertwining operator on 
symmetric polynomials from the TPDs of Dyson's model and symmetric Dunkl processes. This is 
the result of Theorem 2, Section IV A. In order to gain a better understanding of this result, we 
will make observations on the particular case of quadratic symmetric polynomials in Section IV B. 
We will find that in the A; ^ oo limit, turns all symmetric functions into a function of the sum 
of their variables. We will prove this fact in Theorem 3, Section IV C. Finally, we will make use of 
the behavior of Vk->oo to investigate the limit /c = /3/2— T-ooof the TPD of the symmetric Dunkl 
process (or its equivalent Dyson's model) in Theorem 4, Section V. We will find that in this limit, 
this process becomes deterministic and its final configuration is proportional to the roots of the 
Hermite polynomials multiplied by the square root of the process time. We will also find that the 
form of Vfe_^oo is responsible for making the symmetric Dunkl process independent of its initial 
condition in the zero-temperature limit. 



II. DYSON'S BROWNIAN MOTION MODEL AND DUNKL PROCESSES 
A. Dyson's Brownian Motion Model 

Dyson's Brownian motion model [1] is a stochastic process in which N particles undergo Brow- 
nian motion in one dimension while interacting repulsively with each other with a strength pro- 
portional to the inverse of the distances between them. Let us denote the Laplacian operator 

by 

A<"-Ei^- (3) 
Dyson's model is described by the Kolmogorov backward equation (see, e.g.. Remark 1 in Ref. [5]) 

-P,(f,y|a.) = -AWp,(t,y|a=) + f E E -^^-d^M^M^)^ (4) 

i=l j=l: * J * 
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where Pp{t,y\x) is the TPD that the particles reach the positions y = (yi, . . . ,yjv) after a time 
t, given that they started from x = {xi, . . . ,xn)- While we denote vectors in using boldface 
symbols, we will use the notation = x ■ x to denote their squared norm. 

The eigenvalues C = {Ci}i=i,..;N of the Gaussian ensembles of random matrices obey the joint 
probability density function (Theorem 3.3.1 of Ref. [3]) 

Pr:ohNAO expi-pn) (5) 

with 

^ = 1^'- E iogio-Cii- (6) 

l<i<j<N 

Because Prob7v,/3(a^) has the form of the partition function of a log-gas in a harmonic potential in one 
dimension after integrating out the kinetic part, one can consider /3 as a real positive parameter by 
giving it the meaning of inverse temperature in this log-gas model. In the same way, one can treat 
/3 as a positive real number in Eq. (4) by regarding /3/2 as the coupling constant of the interaction 
(drift) term. Therefore, when P is extended from the three values 1,2 and 4, it carries two different 
physical meanings. Loosely speaking, the inverse-temperature interpretation of P still holds in 
Dyson's model: when /3 — )• 0, the interaction term in Eq. (4) vanishes leaving an A?^-dimensional 
Brownian motion, which means that the fluctuations in the system are much larger in magnitude 
than the interaction between particles. This is the analog of a high temperature condition. On the 
other hand, when /? — t- oo, the diffusion term is negligible compared to the interaction term, and 
so the randomness in the system disappears, which is the analog of a zero-temperature condition. 

A consequence of the results in Refs. [43] and [36] (see, e.g., Ref. [38]) is that, for A^-dimensional 
vectors x and y ordered so that yi < yj and Xi < xj for i < j, the TPD of Dyson's model is given 

by 



Nie-{xHy')/2t ^ r r(l + /3/2) T / \ ^3 



p^{t,y\x)= -^^^^^^^^ n 



r(i + W2) 



(2//3) X y 



where T{x) is the gamma function, h]sf{y) = Y\i<i<:j<N{yj ~ Vi) the Vandermonde determinant 
and oJ^Q^^^^ is the generalized hypergeometric function (see Appendix A). In the case /S = 2, the 
non-colliding Brownian motion, Dyson's model has the following TPD as shown by Grabiner [8]: 

PB-2it,y\x) = det {- = . (8) 

^-^^'^"^ hN{x) i<i,j<N [ / 
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B. Dunkl Processes 

Within our present setting, we are interested in the ^-type Dunkl operators given by Eq. (1) 
(see Appendices B and C). Using these operators, one can define Dunkl processes as follows [36]: 
consider the TPD, pk{t,y\x), that solves the Dunkl heat equation 

d 1 ^ 

-Pk{t,y\x) = -J2Thk{t,y\x) (9) 

1=1 

as its Kolmogorov backward equation. Then the stochastic process that obeys pk{t,y\x) as its 
TPD is defined as the Dunkl process of type A. Direct insertion of Eq. (1) into Eq. (9) yields 

d 1 ^ ^ k d 

—Pk{t, y\x) = -AWpfe(i, y\x) + ^~^.d^M*^ 2/1*) 

1=1 j = l: * .? ' 

fe Pkjt, y\x) - Pkjt, vlcTjix) 

i=l j=l: ^ 1 ^' 

Reading ofi' the terms of this Kolmogorov backward equation gives us information about Dunkl 
processes: the first term, a Laplacian acting on cc, is a diffusion term, while the second term is a 
drift term that drives the process away from the surfaces defined by the equations xi = xj for any 
i ^ j. The third term is a spontaneous jump term, which instantaneously makes the process jump 
from X to some permutation of itself so long as Pk{t, y\x) is not symmetric under permutations of 

X. 

With the observation about the jump term in mind, the comparison of Eqs. (4) and (10) reveals 
a clear similarity between the two Kolmogorov backward equations. Let us define the symmetric 
group Spf as the set of all the permutations of N objects with composition as its group operation. 
The permutation p G Sn of the vector x is given by 

= (xp(i),...,a;p(7v)). (11) 

We will consider the TPD of a Dunkl process with a symmetric initial condition. We will refer to 
such a process as a symmetric Dunkl process. Given the symmetric distribution 

f^Uz) = Yl '^'^(^ - P"")' (12) 

where we define S^{x — y) = YljLi ^{xj — yj), we write 

Pki't^y\^)= Pk{t,y\x)i^Uz)d^z = y2 Pk{t,y\px), (i3) 
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with d^z = d^^j- We show that the Kolmogorov backward equation obeyed by the TPD of 



the symmetric Dunkl process, p^{t,y\x), is identical to that of the TPD of Dyson's model in the 
following theorem. 

Theorem 1. The TPD of a symmetric Dunkl process, solves the Kolmogorov backward equation 
of Dyson's model with parameter 



13 = 2k. 



(14) 



That is, the equivalence 



(15) 



is established. 

Proof. Comparing Eqs. (10) and (4), we see that the only difference between the two is the third 
term on the RHS of Eq. (10). The bijective nature of permutations allows us to write 



52 



^=1 Kii) ' 



N N Q N N 

EE = EE 



J= 



and 



N N ^ 

EE ' 



N N 

EE 



1 ~ ^p(i)pU) 

1 j=i ~ ^-j)^ i=i j=i (^pw ~ ^pu)y 



for any permutation p G S^. Therefore, 



^ ^ 5 k^^ l-a^ 



1=1 j=i ^ i=i 7=1 



3= 



E Pk{t,y\px) 

peSN 



= E 



lA 52 



k 



O k ^ - ^ - 1 — (T 



EE 



N N 



jOfe(i,y|pa;) 



9 d 



(16) 



and in particular, since p|(t, y|a5) is symmetric in x, the third term in the first line vanishes, leaving 
us with Eq. (4) provided we set k = P/2. □ 
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Note that, because Pk{t, y\x) is a probability density function on y, its integral with respect to 
y over is equal to one, and so p%{t, y\x) integrates to N\. A consequence of this theorem is that, 
because both and p| follow the same Kolmogrov backward equation, they must be equivalent. 
We will exploit this fact in Section IV. 



III. THE INTERTWINING OPERATOR 



In the previous section we obtained some information about the TPD of Dunkl processes and 
how it is related to Dyson's model. However, this TPD can also be calculated using the intertwining 
operator T4, Eq. (2) [31, 42] (see Appendix D for details). While we will treat Vk from a mostly 
general perspective in this section, we will assume that Vj- is the A-type intertwining operator. 

Dunkl originally defined Vk for the case in which the function f{x) in Eq. (2) is a polynomial, 
and proved its existence for A; > using a recursive construction. However, its explicit form is 
unknown in general, and known only in a few particular cases. Prom Eq. (1), one can note that 
when k = 0, the intertwining operator is the identity operator. 

The intertwining operator relates A^-dimcnsional Brownian motion to Dunkl processes as follows. 
Let us consider the A'^-dimensional heat equation 

(^-^A(-^)po(i,l/W = (17) 



with its solution 



given the initial condition pQ{0,y\x) = 5'^{x — y). This solution is both the Green function of the 
diffusion equation and the TPD of the AT-dimensional Brownian motion. When we apply Vk on 
this equation, we obtain 

^^^^ll-^^^'^O^"^*'^''"^^ (^"^^^')^'^°^*'^''^''^ ^^^^ 

because of Eq. (2) and the fact that V^ commutes with the partial derivative relative to time. From 
the above, we deduce that the Green function for the Dunkl heat equation, which we denote by 
Vk{t,y\x), is given by 

Tk{t,y\x) = VkPo{t,y\x). (20) 

On the other hand, the Green function for the Dunkl heat equation can be calculated directly 
using the Dunkl transform [35] , a generalization of the Fourier transform, in the same way that the 
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Fourier transform can be used to derive Eq. (18), see Ref. [36]. In order to define the Dunkl trans- 
form, one needs to calculate the eigenfunction 0^ of the Dunkl operators {rj}j=i^..._jv of parameter 
k with eigenvalue ^, 

Ti(f>k{x,^) = ^iMx,^), i = l,2,...,N, (21) 

subject to (t)k{0, ^) = 1 with = (0, . . . , 0). When k = 0, Eq. (21) is 

d 

— (/>o(a;,4) = ^ifpoix,^), 

and its solution is 

Mx,^) = e^--. (22) 

This eigenfunction is used as the kernel of the A^-dimensional Fourier transform with x replaced 
by ^/^x. With Eq. (22) in mind, we adopt the form 

Mx,^) = Vke^""- (23) 
Insertion of Eq. (23) into Eq. (21) yields 

TiMx, = TiVke^-^ = Vk^e^-^ = Vk^ie^'^ = C^Ffce^''" = ^i^kix, i), (24) 

where the second equality follows from the definition of the intertwining operator and the fourth 
equality follows from the fact that is linear. In addition, when i^fe(x,^) is written as a power 
series, each term is a homogeneous polynomial of x, whose degree is conserved by V^, and when 
the limit a; ^ is taken, all homogeneous polynomials of degree larger than zero vanish; the only 
remaining term (the zero-order polynomial) is equal to one. Therefore, (t>k{x, as given in Eq. (23) 
is the solution of Eq. (21) with the condition ^^(0,^) = 1. 

In analogy to Eq. (22), (/)fc(a;,^) has the following properties [44, 45]: given a;,t/ G M, 

<t>k{x,y) > 0, 
\(pk{V^x,y)\ < 1, 

(pk{x,y) = (pk{y,x), 
<^fe(Aa;, y) = <t)k{x, Ay), A e C, (25) 

and 



(t>k{px,py) = (t)k{x,y), 



(26) 
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with p G Sn- When the Dunkl operators are not of the A type, (/>fe(a;,y) obeys property (26), 
but the operator p is given by the root system under consideration and Eq. (Bl), and it is not an 
element of Sn in general (see Ref. [37]). Due to these properties, ^fc(£C, ^) can be used as the kernel 
for a transformation similar to the Fourier transform by replacing x by ^J—lx. For this reason, 
hereafter we denote this eigenfunction by 

Ek{x,i) = 4>k{x,$,). (27) 

This function is known as the Dunkl kernel [42] . 

Using the Dunkl transform mentioned above, one can calculate Tk(t, y\x) and the TPD of Dunkl 
processes in general (see Appendix D, Eq. (D24) for the explicit relationship between Ffe and Pk)- 
However, because the Dunkl kernel is contained in it, and the intertwining operator is unknown in 
general, this TPD is not explicitly specified. In our case, pk{t^y\x) is given by 

The last factor on the right is the only unknown part of this TPD, as mentioned before. Let us 
note that in order to obtain pk. we need the Dunkl kernel, which in turn requires us to calculate 
the effect of the intertwining operator on the exponential ^ . Conversely, if we have pk, we can 
obtain the Dunkl kernel, and by expanding it into sums of suitably chosen polynomials, we can 
extract the effect of the intertwining operator on any polynomial. 

IV. EXTRACTION OF Vk FROM DYSON'S MODEL 
A. Theorem 

Here, we will obtain an expression for the effect of Vk on symmetric polynomials. Using Theo- 
rem 1, we combine Eqs. (13) and (28), equate the result to Eq. (7), and rescale x and y by a factor 
of \/t to obtain 

Ek {px, y) = Nlo^'^"^ {x, y) . (29) 

P&Sn 

Prom Eq. (29), we extract the effect of the intertwining operator when applied to a symmetric 
polynomial. As preparation, let us introduce the following definitions: consider an arbitrary integer 
partition, /x, as defined in Appendix A. Let represent the multiplicity of the j-th (distinct) part 
of /X, where the subscript P in Zp refers to the number of distinct parts of p,. In the cases where 
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l{fi) < N, there are N—l{fi) zero parts in /x and therefore Ip = AT— /(/x) accounts for the multiphcity 
of zero parts in the first N parts of /x. For example, if AT = 6 and /x = (5, 3, 2, 2) = (5, 3, 2, 2, 0, 0), 
then P = 4 and = 1, I2 = 1, = 2 and Z4 = 2. Using this notation we define the following 
multinomial coefficient: 

Mif,,N)= (30) 
n • ■ ■ ■ 'p- 

This function represents the number of distinct permutations of fi when it is considered as an 
A^-dimensional vector. See Appendix A (in particular, Eqs. (A3), (All), (A16) and (A17)) for the 
definitions of the monomial symmetric functions mx{x), Jack functions 'P^^''\x) and the various 
quantities involved in the theorem below. 

Theorem 2. The action of the intertwining operator on a symmetric polynomial is given by the 
equation 

r:Kr)<N "^r^W 1^) (fcTV)^ ' > 
\r\=\\\ 

That is, the intertwining operator of type A maps the monomial symmetric functions mx{x) onto 
a linear combination of Jack functions of parameter 1/k, Vt^^''\x). 

Proof. We first expand the LHS of Eq. (29) in terms of symmetric polynomials. For this purpose, 
we expand the symmetrized exponential X^pg5^ exp{px ■ y) in terms of the monomial symmetric 
functions m^(a;). 

oo ^ 

E "EE E ^ E U(-Mnr-" 

|iu|=n T(/i) distinct 



it! 



E E (En-rln!/;"" 

lx:l{lj)<N ^' t&Sn- l^pe5jvj=l )j=l 



r(//) distinct 



= E ;^(e nv''"''l E n^r"' p^) 

t(/i) distinct 

In the last line we have used the substitutions j' = p{j) and p'(j') = T[p'^{j')]- The last term on 
the right is, by definition, m^(y). The term inside the braces is equal to mn{x) multiplied by the 
number of non-distinct permutations of ji. Using Eq. (30), we write the above as 
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Applying on the above yields the LHS of Eq. (29) expanded in terms of a sum of monomial 
symmetric functions. The next step is to eliminate the variable y using the orthogonality of Jack 
functions. Insertion of Eq. (A12) in Eq. (33) after applying Vk yields 

Using Eqs. (34) and (A19), Eq. (29) becomes 

IJ,:l{n)<N^ ^ ^ i/:i'</i t:1{t)<N ^ )t 

and from the orthogonality of Jack functions [29] and the linearity of Vk, which acts only on a;, we 
can equate the coefficients of the same Jack functions of y to obtain 

M=\t\ 

This relation is solved for Vkmx{x) if we apply the sum J2t "^rxi^ / k) , which yields 

Vkmxjx) ^ cAl/k)vi'/'\x) 

XmiX,N) - ^^2^^^-^A(V^)^, (1/,) ^^^^(V.) . (37) 

M=1A| 

Prom this last operation, the theorem follows. □ 

Prom Theorem 2, we can recover known results from the cases k = (where Vj- must reduce to 
the identity operator), k = 1 (which should yield the known TPD of the non-colliding Brownian 
motion) and the case in which all particles start from the same position. We rederive these results 
in detail in Appendix E. 



B. Effect of Vk on symmetric quadratic equations 

Although Theorem 2 is consistent with known results, it does not give an immediate under- 
standing of the effect of Vk on symmetric polynomials. The purpose of this subsection is to consider 
a simple case of Theorem 2 and make observations of the general properties of Vk in order to ob- 
tain insights on the behavior of symmetric Dunkl processes and Dyson's model. Let us consider 
the effect of the intertwining operator on the quadratic equations 7/12(0;) = J2^i^i — ^ 
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Till (a;) — 5^i<j<j<iv ^j^j — 1- 



^ A: + l ^ 2k 

l/fc 2^ 2(A;iV + l)^''^+ A;iV + l ^ a;,a;, -1 (38j 

l<i<j<N ^ ' j=l l<i<j<N 

Prom the coefficients of tlie sums above, it is clear that when A; = the polynomials on the LHS 
are unchanged. When k = 1, these equations become 

^ 2 ( ^ \ 

Vk=iJ2x] = J— J Y^x] + Yl ^^^j- = 1' 

j=l ^ \j=l l<i<j<N ) 

Vfc=i E = %^ + n E + jvTi E ^i^i = i- (39) 

l<i<j<Ar ^ > j=l l<i<j<N 

If we take the limit A; — > oo, the two quadratic equations become 

^ 1 ^ 2 l/^V 

^feE^r-^]vE^i + ]v E ^^^^■ = ]v E^^- 
Vk E^j + E = ( E^i I =1- (40) 

l<i<j<N j=l l<i<j<N \j=l 

These equations represent two pairs of planes, namely 

N N 



I 2N 

J2 Xj = ±^^N and ^ Xj = ±^ 



(41) 



i=i j=i 



We plot Eqs. (38) for the particular case AT = 3 in Fig. 1. The line defined by the vector (1, 1, 1) 
is represented by the diagonal line in each of the subfigures. Prom the figure, we can see that as 
the value of k increases, the surfaces get stretched equally in all directions orthogonal to the line. 
In the limit k — t- oo, the surfaces become planes, and there are two of them because we have chosen 
quadratic polynomials, i.e., |A| = 2. 

Prom these observations, we expect that the effect of V^^oo on the monomial symmetric func- 
tions m\ is to reduce them to the |A|-th power of mi (a;) = J2jLi ^j- The power must be |A| because 
Vk conserves the degree of all homogeneous polynomials by definition. We prove this fact in the 
following theorem. 
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(a) \4=oE-.i^.— 1 



(b) V^fe^iE =1^.-1 



(c) Vk^^ E ■=! = 1 




(d) Vfe=o El<i<j<3 ^i^i = 1 (e) Vfc = l El<i<j<3 ^i^i = 1 (f) 'l^fc-^oo El<i<j<3 ^i^j = 1 

FIG. 1: Effect of the intertwining operator on the equations m2{x) = 1 and mii{x) = 1. 

C. The fc — > oo limit of Theorem 2 



Theorem 3. The limit k ^ oo of Theorem 2 is 



Hm Vkmx{x) 

k—^oo 



M{X,N) 



N 



(42) 



Proof. The quantity 

Cr{l/k) 



n 



{Ti-j + k{T'-i + l)) 



d,{l/k){kNf"'^ (^(^ - i + 1) + i - 1)(t, - J + 1 + A:(rj - i)) 

tends to zero as A; — t- oo in generaL The only case in which this quantity does not vanish is 
when rj = i for any (i,j) G r. Now, rj represents the number of parts greater than or equal to 
j in T (see Appendix A), so we can conclude that this condition is satisfied only when rj = 1, 
i.e., T = (ri,0, . . .). In other words, only partitions of length equal to one satisfy this condition. 
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Therefore, 







whenever i(r) > 1 and 



<(l/A;)(A;iV)i'/*^^ (^^^ + i - - J + 1) /iAr(ri-j + l) iV^n! 

when /(r) = 1. Then, Vk becomes 

due to the fact that if Z(r) = 1, then r' = (l,l,...,l), a partition composed of ones with |A| parts 
(see Eqs. (A4) and (A6) for the definitions of a|A|,A ei(a;)). We finish our calculation giving an 
explicit form of the matrix components a|A|,A by expanding (ei(a;))l^l in terms of mr(£c). We write 

V" lAI- 

(ei(a;))l^l=^xj = ^ —mr{x)= ^ a|A|,rm^(a;), 

\j=l J t:1{t)<N ' t:1{t)<N 

\r\=\X\ \t\=\X\ 

and deduce that 

«|A|,A = ^, (43) 
which proves the statement. □ 
Note that the equation 

M{X,N) 



E^j) =1 (44) 



clearly represents one (for |A| odd) or two (for |A| even) planes with normal vector equal to 

1 = (1,...,1). (45) 

This is consistent with Figs. 1(c) and (f). 



V. ZERO-TEMPERATURE LIMIT OF DYSON'S MODEL 



The physical implications of the observations in the previous section involve the behavior of the 
symmetric Dunkl process when k tends to infinity. This limit is equivalent to the zero-temperature 
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limit (/3 — >■ oo) of Dyson's model. Although the observations in the previous section suggest that 
the limit limfe_>.oo X^pgSj^ y) exists, one must be careful when taking this limit in Eq. (7). 
In particular, the fact that k = P/2 is the coupling constant of the interaction between particles in 
Dyson's model implies that this interaction separates the particles infinitely fast when A; ^ oo. In 

order to obtain an intuitive picture of what occurs at zero temperature, we consider the Kolmogorov 
forward equation associated to Eq. (4) with the parameter k in accordance with Theorem 1: 



N N 



pt:it,y\x) 



ipvhji- yjf ' 



(46) 



We will now apply the scaling y s/kv. Under this scaling, Eq. (46) becomes 



§-pUt,Vkv\x) = lA(Vfe(t,Vfct;|a;)-X; 



AT AT 



2k 

I Y^^ Pl{t,Vkv\x) 

i=l j=l VI J J 



1=1 j=l ■> 
3+i 



dvi 



(47) 



Furthermore, if we assume that limfe_^oop|(f, \/fct;|a;)A;-'^/^ = p^(t, v|a;), taking the limit A; — >■ oo 
transforms Eq. (47) into 



N 



^ pl^{t,v\x) = -J2 ^ 



dt 



dvi 



P%o{t^v\x) 



N 

^ Vi- Vj 



(48) 



a Kolmogorov forward equation with a drift term but without a diffusion term. This equation 
suggests that the scaled process on v is actually deterministic. Physically speaking, this is a 
natural consequence of taking the k ^ oo limit, since the randomness of Dyson's model is removed 
at the zero-temperature limit. Therefore, the argument above implies that p%{t, ■\/kv\x)k^^'^ should 
converge to a function of the form 



lim pl{t,Vkv\x)k^/^ = S^[v - F{t,x)], 

fe->oo 



(49) 



where F{t,x) denotes the vector solution of a deterministic equation of motion. 

Interestingly, the evolution of the scaled process considered above is related to the Hermite 
polynomials defined by 



dx^ 



(50) 



17 



More specifically, the process described by Eq. (48) depends on the N roots of Hn{x), here denoted 

by 

zn = {zi,n,---,zn,n), (51) 

and their permutations pzj^, p € S'at, where Zi^M is the i-th root and Zi^jsj < zj^jsf for i < j. These 
roots are known to be real [47]. Also, because the Hermite polynomials obey the relation 

Hn{-x) = {-1)^Hn{x), (52) 

their roots add up to zero, i.e., 

N 

J2hn = 0. (53) 

We address the zero-temperature limit in the following theorem. 

Theorem 4. The — > oo limit of the TPD of the symmetric Dunkl process p^^{t,^/kv\x)k^/'^ is 
independent of x and is given by 

pl,{t,v\x)=pl,it,v) = ^ 6^[v-V2ipzN]. (54) 

peSN 



Proof. We first calculate the limit 



lim V Ek{px,y). (55) 



fe— >oo 

P&Sn 



Using Theorems 2 and 3 and the expansion in Eq. (33), we obtain 



= iV!expf"" ^fa-^'V (56) 



N ^ 



N 

We now turn our attention to pl.{t,y\x). The probability that the symmetric Dunkl process 
reaches an infinitesimal volume ch^ y around y from x after a time t is pp,{t,y\x) d^y. This 
probability, after the scaling y — )• ^/kv equals 

pl{t, y\x) d^y = p^t, \^v\x)k^^^ d^v = exp { log (p|(t, Vkv\x)k^/'^) } d^v. (57) 
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When k is sufficiently large one can use Stirling's approximation [48] combined with Eqs. (13), (28) 
and (56) to approximate the probability density in Eq. (57) as the exponential of 



k 



N 2 



y (iV - 1)(1 - \ogt) - jlog j + 21og \hN{v)\ - 

j=l 



X 



2 



+^{x .l){vl) + -\ogk + log(7V!) - - log(27ri) - -. (58) 

Let us denote the term in the square brackets by Fi^{v,t). According to this expression, in the 
limit A; — >■ oo the probability density goes to infinity whenever Fj^{v,t) is positive, and vanishes 
when it is negative. The lower-order terms in k are negligible in either of these two cases. However, 
this distinction is difficult to make due to the presence of the Vandermonde determinant of v in 
FN{v,t). We show in Appendix F that FN(v,t) is negative except in the points 

V = V2tpzN- (59) 

The consequence of this is that the probability density (57) vanishes at v ^ V^iZN or its permu- 
tations. On the other hand, when v = \/2tz^ or its permutations, kFN{v,t) (see Appendix F) 
and the term proportional to s/k vanish, because the roots of the Hermite polynomials add up to 
zero. That is, 

1=1 

The only fc-dependent term that remains is ^ log(A:), which tends to infinity with k. Therefore, we 
can write 

lim pl{t,Vkv\x)k^/'^d^v oc V S[v{t) - V2ipzN]. (61) 

peSN 

However, we know that p^t, \'%v\x)k^^^ d^v is equal to A^! when integrated over M^, for any 
value of k. Similarly, the integral over of the sum of delta functions on the RHS is also equal 
to AT!. Therefore, the two quantities above are equal. □ 

Note that the effect of limfe_>.oo X^pgs^ Ek{px,y) on the final result is that it makes the TPD 
p^{t,v\x) independent of x. When FN{v,t) attains its maximum value, the symmetric Dunkl 
kernel is equal to one for all values of x, which leaves A;^/^ as the dominating term in the limit. 
Therefore, the form of the intertwining operator and the properties of the Hermite polynomials are 
ultimately responsible for removing the dependence of p| on x in the zero-temperature limit. 
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VI. CONCLUDING REMARKS 

We established the equivalence between Dyson's model and symmetric Dunkl processes under 
the parameter relation (14) in Theorem 1. This correspondence yields a physical interpretation of 
the abstract parameter k. The relation /3/2 = k implies that one can give to k the interpretation of 
inverse temperature given to the parameter (3. We have also extracted the effect of the intertwining 
operator on symmetric polynomials in Theorem 2 using Theorem 1. It is of interest to note that 
the functions involved, which are commonly found in combinatorics and representation theory 
[28, 29] and in quantum integrable systems [32-34], are found here through the comparison of two 
multivariate stochastic processes. This result stresses the importance of the role that symmetric 
polynomials play in the setting of stochastic processes. While it is known that Schur functions 
are essential in the study of the non-colliding Brownian motion (/3/2 = k = 1) [7, 9, 26, 27], their 
role in the general k case is played by Jack functions [4, 43], which makes them equally essential 
for the study of Dyson's model. We showed in Theorem 3 that the limit /c — t- oo of Theorem 2 
exists, and that in this limit has a simplified form. This fact allowed us to calculate the zero- 
temperature limit of the symmetrized Dunkl kernel, and in turn obtain the TPD of the symmetric 
Dunkl process, which is the result of Theorem 4. The positions of the Brownian particles are locked 
to a deterministic path, which is a natural consequence of the lack of thermal vibrations in this 
regime. 

Moreover, we pointed out the following two properties of the symmetric Dunkl process in the 
limit /c — 7- oo in Theorem 4. The first of them is that the final configuration of the Brownian 
particles depends on the roots of the Hermite polynomials rather than being an equally-spaced or 
simpler configuration. This is a consequence of the form of the interaction between particles: the 
Vandermonde determinant of v appears in p^{t,-\/kv\x) as a death factor that makes the TPD 
vanish whenever Vi = Vj, i ^ j, pushing the particles away from each other. Furthermore, the 
partial derivatives of the Vandermonde determinant are responsible for the form of the drift term in 
Eqs. (4) and (10). Therefore, a different interaction would lead to a different form for the function 
F/v(i),t), which would in turn have its maxima in points different from zn or its permutations. 

The second is that in the zero-temperature limit, p|(i, Vkv\x)k^/'^ is independent of the initial 
configuration x. Because the process is deterministic in this limit, one expects it to depend on 
the initial conditions from which it has evolved. Contrary to this expectation, the nature of the 
intertwining operator at A; — > oo and the symmetry that the roots of the Hermite polynomials obey 
take away that dependence. 
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In summary, we have succeeded in extracting an expression for the intertwining operator, we 
have observed its effect on symmetric polynomials and we have used it to extract information about 
Dyson's model and its behavior at zero temperature. As mentioned in Sees. I and II A, random 
matrix theory and its stochastic version, Dyson's model, have been thoroughly studied and applied 
in numerous branches of physics and mathematics. Therefore, in view of the present results we 
expect Dunkl processes and Dunkl theory to be applied more broadly in physics in general, and in 
non-equilibrium statistical mechanics in particular. 

There are several aspects of Dunkl processes that require further study. An expression for the 
intertwining operator for non-symmetric polynomials is desired in order to obtain the TPD of non- 
symmetric ^-type Dunkl processes. This, in turn, would open the possibility of studying closely 
the nature of the spontaneous jumps that these processes undergo due to the jump term in their 
Kolmogorov backward equation (the third term on the RHS of Eq. (10)). Finally, the existence of 
explicit forms for the symmetrized Dunkl kernel (the LHS of Eq. (29)) for the B and D-type root 
systems as given in Refs. [43, 49] should allow for other types of Dunkl processes to be applied in 
physical systems where symmetry with respect to those root systems is required. 
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We begin by listing briefly the definitions involving integer partitions that are necessary to 
properly define the multivariate functions used in Sees. IV and V [29]. A partition r of an integer 
n is a non-negative integer vector such that its components are arranged in decreasing order while 
adding up to n. The components of r are called parts, and its length, Z(r), is the number of 
non-zero parts in it. The sum of the parts is called the modulus of the partition and is denoted by 
|t|. In other words, r is an integer partition of n if Tj > Tj for i < j, and 
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Appendix A: Multivariate Special Functions 



/(r) 




(Al) 



1=1 
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For example, r = {ti,T2,. . .) = (4, 4, 2, 0, . . .) is an integer partition of n = 10 and has length 3. 
Additionally, the symbol r! denotes the product n^^^^Tj!. We will also use the natural ordering 
for partitions of the same integer defined by 

A < r <^ Xi + ... + Xj <Ti + ... + Tj'^j = 1,...,N, (A2) 

where A and r arc assumed to be different partitions. 

We say that the ordered pair is in the partition r (i.e., G r) when 1 < i < 1{t) and 
1 < i < "^i- The number rj is the number of parts r; such that r; > j. Finahy, the partition formed 
by the rj is denoted by r', and it is cahed the conjugate partition of r. 

Let us now introduce three families of symmetric polynomials, keeping in mind that we will con- 
sider polynomials of N variables. Given an integer partition r, the monomial symmetric functions 
are given by 

mr{x) = J2ll^7'''^ (A3) 

a j=l 

where the sum is taken over all permutations a such that each monomial YijLi -^J^'^' distinct. 
For an integer < n < N, the elementary symmetric function en(x) is given by 

n 

i<n<-<i„<Ar j=i 

For example, eo{x) = 1, ei{x) = X^^i^j, 62(0;) = J2i<i<j<N ^i^j eN{x) = Hili^i- When 
the subscript of e is a partition, it is given by 

1{t) 

eTix) = Yleri{x). (A5) 
1=1 

It is known that the relationship between the monomial and elementary symmetric functions is 
given by 

(^r'ix) = J2 (^rxinxix), (A6) 

A:A<T 
|A|=|r| 

where a-j-x is an upper triangular matrix with diagonal elements equal to one. Note that the 
subscript on the LHS is the conjugate partition of r. 
The Schur function St{x) is given by the relation 

.^-^ deti<i j<Ar[x^''''^~*] 

Sr{x) = mr{x) + > Krxmxix) = — — ^ . (A7) 

deti<,,,<iv[a:f 1 

|A|=|r| 
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(The second equality is the Jacobi-Trudi formula [28, 29]). The matrix K^x is called the Kostka 
matrix. 

The Jack function of parameter a > 0, C^\x) is defined as the polynomial eigenfunction of 
the operator [50] 

+ - E -^^\4-Hx) = E^,^4-Hx) (A8) 

with eigenvalue 

AT 

^r,a = Y.^A^o - 1 - 2(i - l)/«) + |r|(iV - 1). (A9) 

There is one Jack function for each partition r, which is unique up to a normalization constant. 
We will use the following two normalizations: the C normalization, which is defined so that 

oo ^ 

E C(")(a;) =exp(xi+X2 + ... + xjv) (AlO) 

n=0 ■ T:l{r)<N 
|r|=n 

for all a, and the P normalization, which is defined by the following linear combination of the 
monomial symmetric functions mx{x), 

r!r\x) = mr{x) + urx{a)mx{x). (All) 

A:A<T 
|A| = |r| 

The partition-indexed matrix Urx{a), as shown in Ref. [29], is an upper triangular matrix whose 
diagonal entries are equal to one, and its non-diagonal entries are sums of ratios of the form 
{aa + b)/{ca + d), where a, b, c and d are non- negative integers. Therefore, it can be inverted, 
yielding the inverse relation 

mr{y)= {u'\x{a)vt\y). (A12) 

\:\<T 
|A| = |r| 

Also, the limits Yun.a^QUTx{pi) and lima-^-oo ""tA (o;) are well-defined. In fact, it is known that [29] 

Vi^\x) = Sr{x) 

V^°^\x) = mr{x), (A13) 

or, equivalently, 

lim Urxia) = a^-A 

a-)-0 

^^ta(1) = KrX 

lim Urx{a) = 5rx- (A14) 

a— >oo 
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Using Jack functions, one can define the generalized hypergeometric function o^o'*^ {x^v) as 
follows [4, 43], 

n=0 ■ T:liT)<N 
\T\=n 

with 1 as given in Eq. (45). For the general definition of the generalized hypergeometric function 
pJ^q°'\ see Ref. [43]. It is desirable to express the above in terms of the P-normalized Jack func- 
tions so that we can apply Eqs. (A13) and (A14). For this purpose, we require the generalized 
Pochhammer symbol {a)^\ which is defined as the product 

^^)r -ii r{a-{i-l)/a) ' 



as well as the functions 



(«J)Gt 

c;(a) = n i^in - i + 1) + r^j - i). (A17) 



With these definitions, one can write [51] 



Cr{a) 

Insertion of the above in Eq. (A15) yields, with a = 1/k, 



7'(")(1) = ^5(:^;^. (A18) 



oJ-0 {x,y)-X^ ^ —-^^ . (A19) 

n=OT:l{T)<N ' ' {Kly )t 

\T\=n 

Appendix B: Root Systems 

The definition of a root system depends on the reflection operator [37], which is given by 

X ' cx 

= X — 2 a, (Bl) 

a - a 

where a,x e M^^, and a^x denotes the vector obtained by reflecting x across the plane normal to 
the vector a. 
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A root system is defined as a finite set of vectors R which is invariant under reflections along 
its own elements, called roots. That is, G R for all ot,^ e R. R has the property that, for any 
one of its elements ^, E R implies 6 = ±1. It is assumed that none of the roots is a zero vector. 

For every root system, a base can be chosen such that all roots are a linear combination of 
the base vectors with either all coefficients negative or all coefficients positive. The elements of 
such a base are called simple roots, and they divide the root system into the positive and negative 
subsystems, here denoted as and The choice of simple roots is not unique, but the positive 
and negative subsystems they define always divide the root system into two disjoint parts of the 
same size (cardinality). 

The multiplicity function is a function of a root in R that returns a real value and that is 
invariant under reflections along elements of R. That is, 

k{a) = k{ua'OL) ^a, a' G R. (B2) 

The multiplicity function can be understood as a set of parameters (or multiplicities) which are 
assigned to each subset of R formed by roots that are related by a reflection (or composition of 
reflections) along some other root(s). 

It is a well-known fact [37] that the multiplicity function reduces to a single parameter for the 
An^i root system. For simplicity, we will denote this root system by A. Let us prove this fact as 
follows. A is given by 

A = {dij = {ei - ej)/V2 : i,j = 1,. . . ,N;i ^ j}. (B3) 

We will choose the positive subsystem 

A+ = {aij = {ei - ej)/V2 : 1 < j < i < AT}, (B4) 

where denotes the i-th unit base vector. This positive subsystem is generated by the simple 
roots 

ai+i,i = (ei+i - ei)/V2, i = l,...,N-l. (B5) 

Hereafter, we will write ctq..^. = o"jj. We must note that the effect of aij on an arbitrary vector x 
is that of exchanging its i-th and j-th components. To see this, we will compute the Z-th component 
of aax: 



{(7ijX)l =Xl- {Xi - Xj){Sil - 6jl). 



(B6) 



25 



It is easy to see that xi remains unchanged for I ^ and that {aijx)i = Xj and {aijx)j = Xj. 
Therefore, the group generated by the reflections along the elements of A with composition as the 
group operation is the symmetric group Sn- 

In view of this property of A, we sec that we can obtain any root from at most two reflections 
of any other root, as shown below. Consider an arbitrary root ctij and apply to it the reflection 
amj , with m arbitrary. This reflection exchanges the j-th and the m-th components of ctij , leaving 
us with OLifYi- If we reflect once more using (jj;, with I arbitrary, we obtain olijyi as desired. Since k 
is invariant under any of these reflections, one obtains 

k{oLij) = k{amj(JiiOiij) = k{cxim) 

in general, and therefore it can be concluded that k is independent of its argument, so it is a single 
parameter. 



Appendix C: Dunkl Operators and Dunkl Heat Equation 

Dunkl operators are given by [30], 

T,f{x) = d,f{x) + Y: fc(a)^^^^^-4^« • (CI) 
oieR+ 

The first term is a derivative along the vector ^. Each of the summands in the second term is 
proportional to the odd part of / along the vector ot. 

Dunkl defined this operator in order to study multivariate orthogonal polynomials and special 
functions related to reflection groups. Once a root system and multiplicity function are chosen, the 
operator commutes with T^i and therefore has many of the characteristics of partial derivatives. 
Hereafter we will use the notation Tj = Tg. as in Section II. 

In Ref. [36], Rosier found the TPD pk{t,y\x) that solves Eq. (9) as its Kolmogorov backward 
equation, with Tj given by Eq. (CI). The generalized heat equation (9), or Dunkl heat equation, 
is given explicitly by [30] 



)~Pk{ 

{a. ■ x) 



d_ 



-l^Koc)^ 7-—r2 • (C2) 
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Appendix D: Dunkl Transform and Transition Probability Density 

The definition of the Dunkl transform requires both the Dunkl kernel defined in Section III and 
the weight function 

Wkix) = Yl Ict-ajl'^^"). (Dl) 

a&R 

For fixed R and k and functions / € L^{W^ ,Wk), the set of integrable functions with respect to 
the weight function Wk{x), the Dunkl transform is defined by the equation [35] 

fm = -! f{x)Ek{-i^,x)wk{x)d''x, (D2) 
Cfe Jrn 

and 

Ck= [ e-«'/2u;,(0d^e (D3) 

is a normalization constant obtained from a Selberg integral [3] . Due to the properties of the Dunkl 
kernel it is easy to see that 

rf/fe(0 = -/ [-ixjf{x)]Ek{-itx)wk{x)d''x. (D4) 
On the other hand, since 

■^Wk{x) = 2 aj'^^Wkix), (D5) 



cxeR+ 



a direct calculation shows that 

/ [Tjf{xMx)wk{x)d^x = - [ f{x)[Tjg{x)]wk{x)d^x (D6) 
for / and g integrable with respect to Wk- Using Eq. (D6), one can prove directly that 

-[ [T^f{x)]Ek{-i^,x)wk{x)d''x = i^jm). (D7) 



The inverse Dunkl transform is defined as 



f{x) = -f m)E,{itx)wkii)d''^, (D8) 

Ck JRN 

and it is known that the inverse transform of the transform of / is equal to / almost everywhere. 
In view of the identity 

^/_^/(«e«^e-.^d"e = /(x + j,), (D9) 
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where / is the A?^-dimensional Fourier transform of /, the generahzed Dunkl translation Ty is defined 
as 

ryf{x) = -f fk{^)Ek{ix,^)Ek{iy,^)wk{^)d''^, (DIO) 

which reduces to Eq. (D9) when A; = 0. Because Eq. (DIO) does not change when x and y are 
exchanged, we see that Tyf{x) = Txf{y), and because the form of Eq. (DIO) resembles that of 
Eq. (D8), we can write 

- / Tyf{x)Ek{-ii, x)wk{x) d^x = m)Ek{iy, ^). (Dll) 

Cfc JRN 

Finally, we list a useful consequence of Theorem 3.10 in Ref. [42] which is valid for y, z G C-^: 

- [ Ek{x, y)Ek{x, z)e-^'/2^fc(x) d^x = e^y"+'"^/^ Ek{y, z). (D12) 

Now we will use the Dunkl kernel, transform and their properties to find the Green function F^ 
of the Dunkl heat equation (C2). First, we find a non-trivial solution to Eq. (C2) in L?-(R.^ ^wj^), 
which we denote gk{x,t). Applying the Dunkl transform on Eq. (C2) and using Eq. (D7), we have 

d 1 ^ £2 

g-^gkitt) = -^J2^i9k{x,t) = -^gk{x,t), (D13) 

i=l 

where gki^,t) is the Dunkl transform of gk{x,t), and integrating with respect to t we obtain 

5fc(^,^) = Ae-*«'/^ (D14) 

where A is a constant independent of t. Let us define the sum of k{ot) over the positive subsystem 
as 

Using Eq. (D12), one finds the inverse transform of this function to be 

g,{x,t) = - [ Ae-^'/'E,{i^,x)wk{^)d''^ = ^^^e-^V{2t). (d16) 

A direct calculation shows that 

\f2T^gk{x,t) = '-^(^-N-2^=l-^gkix,t). (D17) 
j=i \ ' 

One could think that the calculation of the Green function is finished by giving g^ a displacement 

y to obtain 

gk{y-x^)^ [ e-*^V2,-^ex,.eyd^^. 
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However, this function clearly fails to solve Eq. (C2). Instead, we must use the generalized Dunkl 
translation Eq. (DIO), 

r-^9kiy,t) = - f Ae-'i"/^Eki-ix,^)Ek{iy,^)wki^)d''^, (D18) 
Cfc Jrn 

which automatically solves Eq. (C2) on both x and y. With the substitution ^' = y/t^ and using 
Eq. (D12) we obtain the Green function 



X y 



m,y\x) = T.^gk{y,t) = [vt'TtJ' ^^^^^ 

In other words, the Dunkl transform of a Gaussian centered at the origin is a Gaussian, but in 
order to displace the Gaussian correctly and obtain the Green function above, one requires the 
generalized translation in Eq. (DIO). 

Note that, since we are interested in calculating a TPD based on F^, it must be normalized. In 
particular, we must impose the condition that, for an arbitrary starting point x and time t, 



I pfe(t,y|a;)d^y = l. (D20) 



As opposed to the case of the AT-dimensional Brownian motion, this TPD is not given by F^. For 
example, if we try to normalize directly, we obtain for the case x = 

1=1 m,y\0)d^y = -j^ [ e-yy^'d^y = ^^^. (D21) 

This leads to A = P /{2Tr)'^/'^ , which is a clear contradiction because of the presence of t. Therefore, 
we realize that the normalization of Tj^ must be calculated with respect to in order to obtain a 
consistent result for A: 

1 = J^^ Tk{t, y\x)wk{y) d^y = A^^ J^^ e-y'/^'E^ ^) w^iy) d^y = Ack- (D22) 
For this result, we have used Eq. (D12). This yields A = c^^ and 

In view of Eq. (D20) , we conclude that 

Pk{t,y\x) = wkiy)Tk{t,y\x) = wk{y)VkPo{t,y\x) = wk ( ^ 1 — '^i^t^ — Vv^' V^/ ' ^^^^^ 

The factor Wk {ij/\/t) can be interpreted as follows: in Eq. (C2), the denominators in the second 
and third terms on the RHS denote a repulsion from the planes defined by ck • a; = 0, a G i?, and 
therefore, the probability of any process finishing at any point such that a • y = should be zero. 
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Indeed, W};{y) = Hagii 1^ ' yl'^^"^ cancels whenever o: • y = for any a E R, which reflects this 
fact. 

The existence of Vk is powerful enough to yield all of the above results with relative ease. 
However, the fact that the explicit form of this operator is unknown imposes a limit to its practical 

applications. 

For the purposes of this work, we are interested in pk{t,y\x) for the root system A. In this 
case, the Dunkl operator in the direction ^ is given by 

T^fix) = d,f{x) + k Y: ^-^^^^PP^ii^-irl (D25) 

l<i<j<N -I ' 

from which Eq. (1) follows. Let us calculate the quantities that form in this particular case. 
First, the weight function becomes proportional to the Vandermonde determinant h]\f{x), 



l<j<l<N 



■^{ei-ej)-x 



(D26) 



2kN{N-l)/2 ■ 



The normalization constant Ch is then 

N 

2kN{N-i)/2 J^^- I'-jYv-zi - - - 2fcJV(JV-i)/2 n r(l + A;) ' 

where we have used a particular case of the Selberg integral (equation (17.6.7) in Ref. [3]). Inserting 
the above in Eq. (D24) yields Eq. (28). 



. 1 / ,-V2|. (^M^k.N^_ (2vr)^/^ fr r(l+jfc) 



Appendix E: Reproduction of Known Results 
1. The Intertwining Operator for A; = 

From Eqs. (1) and (2), one sees that should reduce to the identity mapping when A; = in 
Theorem 2. We prove this in the following proposition. 

Proposition 5. In the limit A; — >■ 0, the intertwining operator in Theorem 2 becomes the identity 
operator, i.e., 

lim Vkmx{x) = mx{x). (El) 

fc-s-O 

Proof. Consider the quantity Cr(l/A;)/[c;(l/A;)(A;iV)i^/''^], with {kNfr'^^ given by Eq. (A16) with 
a = kN and a = 1/k. 

c;(l/A;)(A;iV)i^/') ^^.f-^^ {k{N - i + 1) + j - l){n - j + 1 + kir'^ - i))' ^ ' 
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The second factor in the denominator becomes 

1 1 '^"^ 1 1 

n r.-j + l + Mrj-i) ^ r,-j + l =n^ = 7! (^3) 

when k tends to zero. The rest of the expression becomes, after reversing the order of the terms 
in the numerator, 

SfT r.-j+A:(ri-^; + l) _'-[^^ J-l + fc(r;_,+l-^ + l) 

^^^^ k(N + +j -1 ^^^^ k(N-i + l)+j-l ' 
1=1 J =1 ^ ' 1=1 j=i ^ ' 

and after taking the zero hmit all terms cancel except for those in which j = 1: 

fc^o-LJ-J-J- fc(iV-i + l) + 7-l iiiv-i + l 

%=i j=i ^ I ^ j—]^ 

The definition of rj is (see Appendix A) 

r'j = #{m\Tm > j} 

for any partition r (here, # denotes the cardinality of a set), and we are now concerned with 

Ki = #{"ikm > n}. 

By the definition of integer partitions, we know that t„i > ti whenever m < i, so t!^_ returns the 
maximum value of m such that Tm = Ti, which we will call /. Here we consider two cases: in the 
first, i ^ I, but then by the definition of partition, i must be less than I and = Tj+i. In the 
second, i = I; then, we know immediately that Tj > Ti+i. Therefore, for each distinct part in r, 
r|. — i + 1 represents the number of parts Tm such that m> i and Tm = Ti. We then conclude that 
the product 111=1 ("^4 — ^ + 1) is the product of the factorials of the multiplicities of the non-zero 
parts of r. On the other hand, a simple calculation shows that 

TT 1 {N-1{t))\ 
1=1 

which has the factorial of the multiplicity of the zero parts of r in the numerator, so we finally 
arrive at 

^^Oc',(lA)(A:iV)(^/'=) t\M{t,N)- 
Using the above and the third of Eqs. (A14) the result follows. □ 
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2. The Intertwining Operator for k = 1 

This case is of interest due to the fact that Jack functions reduce to Schur functions [29], and 
because of the shghtly simphfied form that Vk assumes, aside from the fact that it represents the 
non-coUiding Brownian motion = 2). 

Proposition 6. When k = 1, the symmetrized Dunkl kernel is given by 

E «)=n(i!) E /-'^'•f (E5) 

peSjv i=l t:1{t)<N i-i-i=l^^^ l^^). 

Furthermore, the TPD of the symmetric Dunkl process is equal to the TPD of the non- colliding 
Brownian motion [8], Eq. (8). 

Proof. At A; = 1, Urxil/k) becomes the Kostka matrix K^-x, and the functions Cril/k) and c'^{l/k) 
become equivalent and cancel in Eq. (31). Also, the generalized Pochhammer symbol, Eq. (A16), 
becomes 

(fciV)f/'.l,.,=n(^-^, (E7) 



which transforms Eq. (31) into [52] 



N 



Vkmx{x)\k=i = XlM{X,N) J] U (n^_ ■ f A^ ^^xSrix). (E8) 

T:l{T)<Ni=l ^ I' + Ti)- 

\r\ = \\\ 

Now, from Eq. (33) and Eq. (E8), we obtain 

E,=,{px,y) = ^'-^^(y) E UwzTrL^^rM^) 

peSN ti-l{ti)<N t:1{t)<N 1=1^ 

\r\=M 

T:l{r)<Ni=l ^ M:^(m)<^ 

I/^I=|t| 

T:l(T)<Ni=\ ^ 

The coefficient of the product of Schur functions can be rewritten as 
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from which we obtain Eq. (E5) . As shown in Appendix A of Ref . [9] , the RHS of Eq. (E5) is the 
Schur function expansion of the ratio of determinants in Eq. (E6). Therefore, inserting Eq. (E6) 
in the combination of Eqs. (13) and (28) evaluated at A; = 1 completes the Karlin-McGregor 
determinant [53] for the one-dimensional heat kernel and yields Eq. (8). □ 

3. Dyson's Model: Pairticles Stcirting Prom the Same Position 

In view of the ^-type Dunkl operator, Eq. (D25), we see that when it is calculated along the 
vector ^ = 1, the difference term cancels and 

Tif{x) = difix) (Ell) 

for an arbitrary differentiable function /. Prom this, one can deduce that the intertwining operator 
should have no effect on / along the line defined by 1. Moreover, the Dunkl kernel should reduce 
to a regular exponential in this case. The consequence of this fact on Theorem 2 is the following. 

Proposition 7. The intertwining operator has no effect on the monomial symmetric function 
mx{x) when it is evaluated at x = xqI, where xq G M. Moreover, the TPD of the symmetric Dunkl 
process whose initial configuration is x = xqI is given by 

piit, y\xol) = pUt, y - xom = n T{ltjk) l^^^^^l'' e-("i)V2*. (E12) 

Proof. First, we prove that Vkmx{x) is independent of k whenever £c is a multiple of 1. For this 
purpose, we derive an expression which will prove useful for our calculations. By definition, the 
C-normalized Jack functions satisfy Eq. (AlO), which is equivalent to 

J2 OfHx) = {xi + X2 + . . . + xnT (E13) 

t:1{t)<N 
\T\=n 

for any non-negative integer n. The RHS can be rewritten in terms of monomial symmetric 
functions, 

{X1 + X2 + ... + xnT = Yl ^■mAx)^ (E14) 

v:l{u)<N 
\v\=n 

but at the same time, we can rewrite the LHS of Eq. (E13) using Eqs. (A18) and (All): 

t:1{t)<N t:1[t)<N '^^ ' ' t:1{t)<N '^^ ' ' ii:l{\x)<N 

\T\=n l''"l=" ImH'* 
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Since the monomial symmetric functions are linearly independent, we can equate the coefficients 
of nil, in Eqs. (E14) and (E15) to obtain 

t:1{t)<N ' ' 

\T\=n 

Now we consider Theorem 2 at £C = xqI. 

..™.(M)^A!M(A.«) E ^(f)" 

t:1{t)<N ' ' 

\A=\A 



ir| = |A| 

We insert Eq. (E16) with n = |A| and v = \io obtain 

FfemA(xol) = M(A, 7V)xl-^l . (E18) 
As expected, this quantity is independent of k. Note that this is the same value as 

mA(xol)=4''' E l = M(A,A^)xt'l, (E19) 

distinct 
permutations of A 

which shows that has no effect on m\ along 1. Therefore, 

^ E^{p{x^\),y) = N\e-'^-y. (E20) 

Inserting the above in Eq. (13) combined with Eq. (28) evaluated at a; = xqI completes the 
proof. □ 

This is consistent with previous results [4, 7]. In particular, when xq = and t = l/{2k), the 
TPD in Eq. (E12) becomes 

p|(l/(2*),^|0) =iVl?^ ni^|*,(^)|-e-»'. (E21) 

This is the joint probability density function of the eigenvalues y of a matrix drawn from the 
/3 = 1, 2,4 ensembles of Gaussian random matrices multiplied by a factor of N\, with /3 = 2k [3]. 



Appendix F: Externa and Maximum value of Ff^{v,t) 



Here, we prove two lemmas necessary to complete the proof of Theorem 4. The first concerns 
the location of the maxima of Fn{v, t). 
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Lemma 8. The maxima of the function 



N 



N 



FN{v,t) = ^(iV - 1)(1 - logo - Y.jlogj + 2log\hN{v)\ - ^ 



(Fl) 



are located atv = v2izjv or any of its permutations, where zn is the vector of roots of the Hermite 
polynomial Hn{x), given by Eqs. (51) and (50) respectively. 



Proof. We first calculate the first-order partial derivatives of F/v(v,t) relative to v and equate 
them to zero. 

N 

(F2) 



FN{v,t)= J2 



dv. 



- — = 



Vi - V-i 



Hence, the extrema of Fjv(v,i) must obey the relation 

^ 2t 

- Vj 



(F3) 



Furthermore, applying the scaling v = yitz we obtain 

AT 



Zi 



E— ■ 

3'-3+i- 



(F4) 



Therefore, it suffices to solve the above equation for z to find the location of the extrema of 
FN{v,t). The second order derivatives of Fjv(v,i) are 



92 



dvjdvi 



2\og\hN{v)\-- 



_d_ 

dvj 



V- 

^-^ 1); 



Vi 



l-.l^i 



Vi-Vl t 



El: 



7 if « = J, 



(F5) 



if i ^ j. 



The matrix formed by the NxN second order derivatives above is negative definite for all vectors v 
with non-repeating components. To show this, we consider an arbitrary real vector u and calculate 

the quadratic form associated to Eq. (F5). 



E 



Ui 



d ^FNiv,t) ^ 
dvjdvi 



N 



(F6) 



l<i,j<N ■> " i=l l<i<j<N 

Here, the equality holds only when all the Uj are equal to zero. This shows that all extrema given 
by Eq. (F3) are maxima. 

Note that, given z, any of its permutations solve Eq. (F4): 

^ 1 

(F7) 



Zp{i) = X] Z 
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for any p E Sn- Now, we prove that the roots of the Hermite polynomial solve Eq. (F4). One 
can write iJjv as follows: 

N 

HN{x) = 2^ll{x-Zj,N). (F8) 
The derivative of the logarithm of Hn is 

N 



Hn{x) ^j,N 

where 



K{x) = -^Hn{x). (FIO) 



We subtract the i-lh term of the sum from both sides to obtain 

N 



H',i.) 1_ ^ ^ 



H.^(x) X - Zi^N ,4^1 X - Zj^N ' 



Taking the limit x — >■ Zi^N yields 



f ' = lim ^544 (F12) 

2:i,Ar - Zj N ^-^^i,N Hn{x) x - Zi N 



At the same time, 

N 



K{zi,N) = lim = 2^ n (^^.^ - ^J.^)- (F13) 

2;-^2i,JV X — Zi N .-^.i . 

3=1 



We will denote the expression in the limit by 

AT 



3-3¥=i 



X - ZiN 



so that 



lim Hi,N{x) = H'j,{zi,N). (F15) 



With this, Eq. (F12) can be rewritten as 



y 1 1 / H'^jx) Hi,N{x) 

Zi,N - Zj N ^^Zi,N Hi N{x) \X- Zi N X - Zi N 

3=1 

1 f H'^{x) Hn{x) 



lim -^^l^lMll^- ""^-^ . (F16) 

x-^Zi,N Hi^N[X) \X - Zi^N (X - Zi^Nr , 
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In the vicinity of Zi^^, Hn{x) is approximately equal to 

1 



H]Sl{x) K, H'j^{zi^N){x - Zi^N) + ^H'l}{zi^]Sl){x - Zi^Nf, (F17) 



where 



Ki^) = -^Hm{x). (F18) 
At the same time, Hermite polynomials satisfy the differential equation [46] 

H'l^{x) = 2xH'^{x) - 2NHn{x), (F19) 
which implies that H'^{zi^f^) = 2zi^NH'j^{zi^N), so 

Hn{x) H'^{zi^N){x - Zj,Ar)[l + Zi^N{x - zi^n)]. (F20) 
Inserting the above in Eq. (F16) yields 

N 



1 1 f H'^ix)-H'^{Zi,N) ^ ^, , A 

- Zj,N ^-^^i,N Hi^N{x) \ X- Zi^N J 

H%{zi^N) 



rr, , ^ Zi,N = Zi,N, (F21) 

J^]\[[Zi,N) 

which completes the proof. □ 

In the second lemma, we find the maximum value of FN(v,t). 
Lemma 9. The maximum value of F^iv^t) is zero. 

Proof. By Lemma 8, the maximum value of Fj<;{v,t) is located at v = \/2izN. This yields 

N ^ 
FN{V2izN,t) = -{N - 1)(1 + log 2) - ^ jlog j + 21og |/ijv(zjv)| - z%, (F22) 

i=i 

which is independent of t. Let us focus first on the term 

N 

4 = E4iV- (F23) 
The derivative property of the Hermite polynomials [46] allows us to write 

H'^ix) = 2NHn-i{x), (F24) 
and using the product representation from Lemma 8 we obtain 

N N N-l 

E n - ^^•.^) = ^ n - ^^-.^-i)- 
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Since the above equation holds for all x, we can expand in powers of x and equate the coefficients 
of the same order. In particular, the coefficients of x^"^ form the equation 

{N-2) Zi,NZj,N = N J2 Zi,N-iZj,N-i, (F26) 

l<j<jr<JV ^<i<j<N-l 

and if we denote the double sum on the LHS by s^r , we obtain 

N 

SN = j;^^N-i. (F27) 

Now, we can show that S]\f = —N(N — l)/4 by induction. Since H2{x) = — 2, its two roots lie 

at X = ±1^1/2 and therefore S2 = —1/2. Then, the formula is valid for N = 2. Assuming it is also 

valid for A/" — 1, we verify that it holds for N: 

AT N {N-l){N-2) N{N-1) 
- iV^'^-^ = 4 = 4 ' 

as desired. Now, since the Hermite polynomials have either odd or even symmetry, the sum of 

their roots equals zero. Therefore, 

, 2 

0=|>>i,jv| =z% + 2sN = z%-2 ^^^~^\ (F29) 



\i=l J 



which finally yields 



Now, let us calculate 



4^^. (F30) 



21og|/i7v(2iv)|. (F31) 

The main quantity we wish to find is the square of the Vandermonde determinant of the roots of the 
AT-th Hermite polynomial. This quantity is known as the discriminant of the Hermite polynomials, 
and we calculate it here following Szego [47]. Using Eq. (F13), we can write 

(_l)NiN-l)/2 (_l)NiN~l)/2 j^N 

= -^^ U^nM = ^ V(^-i) n^^-iK^)- (F32) 

i=l i=l 

The last equality follows from Eq. (F24). Let us focus on the last product: 

N N N-l N-1 N 

n^iv-i(^i,iv) = 2^(^-1) n n (^^.^ - = ^^^^"^"'^ n n(^^-.^-i - ^^^) 

1=1 i=l j=l J=l i=l 

7V-1 

= n HN{zj,N-i). (F33) 

3=1 



38 

Prom the recurrence relation 

HNix) = 2xHn-i{x) - 2{N - l)HN-2ix) (F34) 
we obtain that Hjq{zj^N-i) = —2{N — l)HN-2{zj,N-i), so the product above becomes 

N N-1 

llHN-liz^,N) = [-2{N - H Hn-2{Zj,N-i). (F35) 

i=l j=l 

For AT = 2, the LHS is equal to 

Hl{-l/^/2)Hl{l/^/2) = -2, (F36) 

for which the general formula 

N N-1 

i=i j=i 

holds. Now, assuming that Eq. (F37) holds for A'' — 1, we verify that it holds for A'': 

N N-l 
HHn-iM = [-2(Ar - H HN-2{Zj,N-l) 

i=l j=l 

N-2 N-l 

= [_2(iV - i)]^-i(_2)(^-i)M/2 Yl f = (_2)^(^-i)/2 Yl f. (F38) 

Therefore, 

(_-,\N{N-1)/2istN ^ 1 ^ 

ihN{zN)f =^-^^rW^^ \{H^-M,n) = i^^^w^,\{3' ■ (F39) 

i=l j=l 

The logarithm of the above is 

^ TV 

21og|^^(zjv)| = ^Jlogj - -(AT - l)log2, (F40) 

and thus the maximum value of Fn{v, t) is zero for all its extrema, due to the fact that Eqs. (F23) 
and (F31) do not change when the roots zn are permuted. □ 
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